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1 Abstract. In this review we discuss possible systematic errors inherent in classical ID LTE 
abundance analyses of late-type stars for the light elements (here: H, He, Li, Be and B). The 
advent of realistic 3D hydrodynamical model atmospheres and the availability of non-LTE line 
formation codes place the stellar analyses on a much firmer footing and indeed drastically 
^ ' modify the astrophysical interpretations in many cases, especially at low metallicities. For the 
, T c ff-sensitive hydrogen lines both stellar granulation and non-LTE are likely important but 
the combination of the two has not yet been fully explored. A fortuitous near-cancellation 
of significant but opposite 3D and non-LTE effects leaves the derived 7 Li abundances largely 
i unaffected but new atomic collisional data should be taken into account. We also discuss the 

p<^ impact on 3D non-LTE line formation on the estimated lithium isotopic abundances in halo 
stars in light of recent claims that convective line asymmetries can mimic the presence of 6 Li. 
While Be only have relatively minor non-LTE abundance corrections, B is sensitive even if the 
■ latest calculations imply smaller non-LTE effects than previously thought. 

Keywords, convection, line: formation, radiative transfer, Sun: abundances, Sun: atmosphere, 
stars: abundances, stars: atmospheres, stars: Population II, Galaxy: abundances 
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1 ■ 1. Stellar model atmospheres and spectral line formation 

• Stellar chemical compositions are not observed: to decipher the spectral fingerprints 
^\ \ in terms of abundances requires realistic models for the stellar atmospheres and the line 
ON ' formation processes. Traditionally abundance analyses of late-type stars have been car- 
ried out relying on ID, time-independent, hydrostatic model atmospheres, which treat 
convection with the rudimentary mixing length theory. All of these are dubious approx- 
imations as even a casual glance at the solar atmosphere will immediately reveal. The 
solar atmosphere, as for other late-type stars, is dominated by granulation, which is the 
observational manifestation of convection: an evolving pattern of broad, warm upflows 
in the midst of narrow cool downdrafts. Because of the great temperature sensitivity of 
the opacity, the temperature drops precipitously as the ascending gas nears the optical 
surface before it overturns and is accelerated downwards. The temperature contrast is 
therefore very pronounced in the photospheric layers, amounting to > 1000 K at the 
optical surface for the Sun (even when averaged over surfaces of equal optical depths 
these rms-diffcrcnccs amount to « 400 K). In addition to ignoring such atmospheric in- 
homogeneities, ID model atmospheres can not be expected to have the correct mean 
temperature stratification because of the simplified convection treatment and the neglect 
of convective overshoot. Even small temperature differences can propagate to very large 
changes in the emergent stellar spectrum because of the non-linearities in the radiative 
transfer and the extreme opacity variations. 
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Over the past decade or so, 3D, time-dependent, hydrodynamical model atmospheres 
for a range of stellar parameters ha ve started to be developed and applied to stellar 
abundance work (e.g. I Asplund 2005 , and references therein). Such 3D models solve the 
standard hydrodynamical conservation equations coupled with a simultaneous solution 
of the 3D radiative transfer equation and the refore self-consistently predict the convec- 
tive and radiative energy transport (see e.g. iNordlund et~aTll200a for further details). 
To make the 3D modelling computationally tractable, the radiative transfer is not solved 
for the many thousands of wavelength points as routinely done in ID model atmosphere 
codes. Instead the frequencies are sorted in opacity and more recently also in wave- 
length space for w 4 — 20 bins for which the radiative transfer is solved. The resulting 
total radiative heating/cooling as a function of atmospheric depth is surprisingly well re- 
produced. The 3D models are based on similarly realistic microphysics (equation-of-state 
and continuous/line opacities) as employed in standard ID model atmospheres. Currently 
there are mainly four different codes b eing used to develop 3D models - S TAGGER (e.g. 
Nordlund et al.l l2009h . c o5bold (e.g. iLudwie, et al.ll2009bh . muram (e.g. IVogler et all 
20051 ) and ANTARES (e.g. iMuthsam et al.ll2009l) - although essentially all abundance re- 



lated work to date has been performed within the first two collaborations. 

Being more sophisticated in the modelling does not automatically translate to being 
more realistic. Over the past few years, substantial effort has therefore been dedicated 
to verify the suitability of the 3D models for quantitive stellar spectroscopy using an ar- 
senal of observational diagnostics. Some of the striking successe s are that the 3D mod els 
accurately predict the detailed solar granulation properties (e.g. Nordlund et al.ll2009l) as 



well as spectral line profiles, including their asymmetries and shifts (e.g. lAsplund et al 



2000), which strongly suggests that the 3D modelling captures the essence of the real at- 



mospheric structure and macroscopic gas motions. Another crucial test is the continuum 
center-to-limb variation, which is an excellent probe of the mean temperature stratifica- 
tion in the solar atmosphere. As shown in Fig.[TJ the latest generation of 3 D models com- 



puted with the STAGGER code reproduces the observations extremely well (jPereira et al 



l2010h : the C05BOLD solar model does a similarly good job (H.-G. Ludwig, private commu- 
nication). This is a remarkable achievement. The 3D solar model not only outperforms all 
tested theoretical ID model atmospheres like the Kuru cz, MARCS and PHO E NIX fl avours 
in this respect, it also does noticeably better than the iHolweger fc Mullen (jl974l ) semi- 
empirical model, which was constructed largely to fulfill this observational constraint. 
The 3D solar model also performs very we ll when confronted with other tests, such 
as the spect ral energy distributio n, H lines ( Pereira et al. 2010l ) and spatially resolved 
line profiles (jPereira et al. 2009a. b). In all aspects the most recent 3D solar models are 
clearly highly realistic and can therefore safely be trusted for abundance analysis pur- 
poses ( Asplund et al. 20091 ). 

The success in the solar case gives some credence to the 3D modelling being similarly 
realistic also for other stars for which far fewer indisputable observational tests are avail- 
able. 3D models now exist for a sizable portion of the HR-diagram covering essentiall y 
spectral types A, F, G, K, M and even later fe.g. lCollet et al.ll2007t iLudwig et al.ll2009bl ). 
Both dwarfs, subgiants, giants and even some supergiants have been simulated and for a 
wide range of metallicities. Qualitatively the granulation behaves similarly to the Sun for 
most of these stars but typically the convection becomes more vigorous towards higher 
T e ff and lower logg. Perhaps the most marked difference compared with ID models are 
apparent at low [ Fe/Hl, where 3D mo dels have much cooler temperatures in the opti- 
cally thin regime ([Asplund et al. 1999). This is a natural consequence of the dominance 
of adiabatic heating over radiative heating in absence of a heavy line-blanketing at low 
metallicities. The low temperatures are bound to affect especially minority species, low 
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Figure 1. Comparison of the o bserved continuum center-to- limb variation as a function of 
wavelength (Neckel fe Labs 1 994f ) against the prediction s for different solar model atmospheres. 
The 3D solar model employed bv lAsplund et all (120091) outperfo rms the ID theoretical MARCS 
model and even the ID semi-empirical iHolweger fe Miillerl l) 19741 ) model atmosphere. 



excitation lines and molecules, with 3D LTE abundance corrections often amounting to 
> +0.5 dex fe.g. lAsplund1l2005l) . One should be aware though that the steep temperature 
gradients may significantly enhance the non-LTE effects compared to the ID case. 

With a stellar model atmosphere it is possible to compute the emergent spectrum, 
which normally for late-type stars is done within the LTE approximation and in ID. It 
should not come as a surprise that this approach often leads to severe systematic errors. 
In recent years, efforts towards rectifying these shortcomings have mainly focused on 
ID non-LTE investigations and 3D LTE calculations, b ut the combination of 3D and 
non-LTE is still largely unexplored (e.g. Asplund 20051 and references therein). When 
solving the non-LTE problem one must solve the rate equations for all relevant levels 
involving both radiative and collisional transitions coupled with a simultaneous solution 
of the radiative transfer equation for all necessary frequencies. Needless to say, it rapidly 
becomes far more cumbersome than in LTE with a great deal of additional atomic data 
that must be considered. While the radiative processes like transition probabilities and 
photo-ionization cross-sections are now in reasonably good shape - especially for the non- 
complicated cases like the light elements - the main uncertainty in non-LTE calculations 
often stems from the poorly known collisional excitation and ionization due to electrons 
and H as well as charge transfer reactions. Especially the inelastic H collisions are a major 
proble m in stellar abu ndance analyses. Mostly simple but likely erroneous classical recipes 
like the lDrawinl ([19681 ) formula are applied, often with an additional ad-hoc scaling factor. 
For a few elements, Li being a nota ble case in point, detailed quantum mechanical or 
laboratory me asurements exist (e.g. Belvaev fe Barklem 2003), which reveal that the 
Drawinl ( 19681 ) formula typically overestimates the real cross-sections by several orders of 
magnitude. An alternative approach is to attempt to calibrate the unknown collisional 
rates using various observations, such as obtaining consistent results from different lines 
or using center-to-limb variations. Not surpri singly, such empirical proced ures suggest a 
range of scaling factors, from ^ 0.1 for Mn ( Bergemann fe Gehren 20081) . to w 1 for O 
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( Pereira et al.l [2009a ) or even > 1 for Fe ( Korn et al. 2003). Whether these are realistic 
values, or merely artificial results due to too simplified modelling (e.g. ID) remains to 
be seen. It is also worth remembering that a common scaling factor for all transitions 
within the same species, let alone between different elements, is highly unlikely. In spite 
of the remaining uncertainties introduced by our still poor handle on all atomic data, 
the alternative of relying on the LTE approximation will most often lead to even more 
severe systematic errors. 



2. The light elements 

2.1. Hydrogen 

As the most abundant element and the main provider of continuous opacity, it is not 
possible to derive the stellar H abundance spectroscopically. The H lines are still of great 
use in stellar physics as a probe of the atmospheric conditions. In late-type stars the 
wings of the Balmer lines reflect the effective temperature and the temperature gradient 
near the continuum forming layers with only a relatively small gravity dependence. The 
great temperature sensitivity arises from the high excitation potential of the lower level 
(Xoxc = 10.2 eV) of the Balmer lines. Fort unately, both the Stark broadening and self- 
broadening are now well established (e.g. Barklem et al. 2000t Allard et al. 20081 ): the 
new broadening results in significantly lower T c ff, especially at low [Fe/H]. 

The H lines are however sensitive to the convection treatment, which in ID normally 
is estimated by the mixing length theory (MLT) and its inherent four free parameters. 
Before comparing the best-fitting «mlt and T e g between different studies it is therefore 
important to bear in mind the various flavours of MLT existing in different ID model 
atmospheres. Given the rudimentary nature of MLT to describe the correct convective 
energy transport and overshoot near/in the optical surface, ID-based T c g-scales from H 
lines will always be unce rtain in an absolute sense even if highly accurate relative val- 
ues can be obtained (e.g. Asplund et al. 20061 ). As explained above, 3D hydro dynamical 
models do provide an attractive alternative through their se lf-consistent computat ion of 



the convection without the need to invoke free parameters. iLudwig et al.l (|2009al) have 
investigated the 3D LTE formation of Ha, H/3 and H'y lines for six different 3D mod- 
els, including the Sun and metal-poor dwarfs/subgiants and compared with ID models 
computed with i dentic al microphysics but wit h different MLT imple mentations; see also 
Sbor done et al. I (l2010h for further discussions. ILudwig et al.l (|2009af) found a very com- 



plex dependence on the 3D corrections to the ID-based estimates of T c s with differences 
amounting to ±300 K depending on the line in consideration and the T e g, [Fe/H] and 
qmlt of the ID model. An extra complication arises from the different line shapes in 
3D and ID, which can not be directly translated to a T c g difference without specifying 
exactly which wavelength regions have been used for the comparison. The prospect for 
using pre-tabulated 3D corrections to ID-based T e g estimates is therefore not encour- 
aging. It would be advantageous to directly compare observations with 3D predictions, 
which should be possible in the near future with the advent of grids of 3D models. 

Until recently it has always been argued that the wings of the H Balmer wings are 
formed in LTE. IBarkleml (|2007j) has investigated whether this is indeed correct and found 
that it is not necessarily so. In particular he found that electron collisions are not sufficient 
to establish LTE. In terms of T e s, the non-LTE calculations would imply 100 K higher- 
values than in LTE. This however depends crucially on the still uncertain inelastic H 
collisions - improved atomic physics calculations are clearly needed to resolve this issue. 
It would also be necessary to perform 3D non-LTE H line computations, especially at 
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low metallicity. Suc h calculations have recently been performed in the solar case with 
encouraging results ( Pereira et al. 201Clt ). 



2.2. Helium 

While in hot stars the very highly excited transitions of He I enable a determination of 
the He abundance (see Kaufer, these proceedings), in late-type stars the corresponding 
atomic levels are not sufficiently populated at typical photospheric temperatures. Instead 
the He I lines (e.g. 1083.0 nm) have a chromospheric origin and can neither probe the He 
abundance nor the photospheric temperature structure. 

2.3. Lithium 

Stellar Li abundance are usually determined only from the resonance line at 670.7 nm but 
in exceptional cases also from the weaker subordinate line at 610.4 nm. Contrary to most 
elements, the simplicity of the Li atom has enabled quantum mechanical calculations 
of all atomic data necessary for non-LTE investigations, especiall y of cross-sections fo r 



collisional excitation and charge transfe r with neutral hydrogen (jBarklem et al.l 120031 ) 



The widely used study of lCarlsson et al.l (|1994j) mapped out the non-LTE effects in late- 



type stars and their detailed variation with metallicity, effectiv e temperature, surface 
gravity, and lithium line strength. Recently, iLind et al.l (|2009al ) revisited the ID non- 



LTE line formation with improved collisional data and treatment of line-blocking (see 
Fig EJ- The investigations have shown that two competing effects govern the population 
levels of Li. Over-ionization from the first excited level of Li I is prominent at low effective 
temperatures due to a strong J v — B v excess bluewards of the photo-ionization threshold 
at 350 nm. When the Li line is weak, the balance between UV over-ionization from 
the first excited state and over-recombination to higher excited levels determines the 
population levels of Li i. The resulting non-LTE corrections are typically minor, smaller 
than ±0.1 dex, but increasing up to +0.5 dex in red giants. However, when the line is 
close to saturation, the loss of line photons establishes an efficient recombination and 
de-excitation ladder, which increases the population of low-excited levels of neutral Li. 
Furthermore, the scattering-dominated resonance line source function drops far below 
B v . The non-LTE corrections then change sign and reach —0.5 dex in extreme cases. 



Lind et al.l (|2009af) have made available convenient routines to interpolate the non-LTE 
corrections for a wide range of stellar parameters and Li abundances. 

In LTE the Li I level populations depend on the local conditions, especially the gas 
temperature. The much cooler mean temperature stratifications in 3D models than cor- 
responding ID models at low metallicity translate to a 3D LTE abundance correction of 
~ —0.3 dex for the L il 670.8 nm resonance line in metal-poor dwarfs and subgiants (e.g. 
Asplund et allll999l ). In addition, the presence of atmospheric inhomogeneities typically 



also increase the line strength. One should be aware however that the steeper tempera- 
ture gradients in the 3D models are likely to boost the over-ionization. Spatially resolved 
observations of the solar surface show how the line actually is weaker in the bright, 
granular upflows, despite the drastic drop in temperature that they undergo. This can 
be understood s ince the intense UV radiation field in the granules produces pronounced 
o ver-ionization (iKiselman 1997[) . 

Asplund et all (|2003l ) performed full 3D non-LTE computations for the Sun and two 



metal-poor stars using a 21-level Li model atom. They concluded that the line strength- 
ening due to the cooler temperatures in the upper atmospheric layers on the one hand 
and li ne weakening from inc reased over-ionization on the other hand largely cancel each 
other. iBarklem et al.l (|2003l ) updated these results in light of new collisional data, espe- 
cially charge transfer reactions and reached similar conclusions. The 3D non-LTE cal- 
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Figure 2. The ID non-LTE abundance corrections for a sel ection of stellar parameters (labelled 
with T c ff / log <?/[Fe/H]) from the study of lLind et a l. (2009a). The solid line denotes the case with 
all transitio ns considered w h ile th e dotted line correspond to the case when the charge transfer 
reactions of iBarklem et af] (I2003T) are ignored. Also shown as dashed lines are the non-LTE 
corrections from lCarlsson et al.l (|!994l ). 

culations thus agree with the ID non-LTE results to within 0.05-0.1 dex. In fact, even 
better a greement between the two is seen in the recent study of extremely metal-poor 
stars by Sbordone et al. ( 2010l) with difference of < 0.03 dex, at least compared with 
the particular family of ID model atmospheres they used. ISbordone et"~aT1 (|2010l ) used 
a smaller atom (8 levels) but computed for more 3D models. They also provided a very 
handy functional relation between the Li abundance and the line strength, which allows 
direct computation of the 3 D non-LTE bas e d abu ndance without invoking any particular 
ID model atmosphere. The lSpite fc Spitel (|l982t ) Li-plateau as inferred from ID analy- 
ses is thus unlikely to be ser iously in danger although the exact slope w ith for example 
[Fe/H] may be affected (e.g. lAsplund et aLll2006l ; TSbordone et alj|2010h . which deserves 
to be studied closer with improved 3D non-LTE calculations. 

2.4. 6 Li/ 7 Li 

The minor isotope 6 Li can be detected through the isotopic shift in the Lil 670.8 nm line. 
The distortion of the line profile is very small and therefore necessitates extremely high- 
quality spectra {S/N > 400, A/AA > 100, 000 and mini mal fringing). At solar metallicity, 
possible blending lines must also be carefully evaluated ([israelian et al.ll2003l ). a problem 
which how ever disappears at l ow [Fe/H]. Th e first positive detecti on in a halo star was 
claimed bv lSmith etail (jl993l ) for HP 84937. fAsplund et ail (|2006l ) boosted the number 
of > 2a detections to 10 (including HD 84937) using high-quality VLT/UVES spectra 
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analysed both in ID and 3D LTE. The derived 6 Li abundance at [Fe/H]< —2.5 can not be 
easily accommodated with standard Galactic cosmic ray production, which has spurred 
a number of studies of more speculative production chann els, including; non-standard B ig 

Bang nucleosynthesis with supersymmetric particles (e.g. Ijedamzik fc Pospelov 2009 ). 

Sev eral potential problems linger over the published 6 Li/ 7 Li analvses. lGarcia Perez et al 
(2009) rave suggested that residual fringing in the observed spectra can cause large un- 
certainties. Their Subaru/HDS s pectra are, however, much more inflicted by fringing 



sp_ 

than the VLT/UVES spectra of lAsplund et al.l (|2006l) . who furthermore carefully as 



sessed their possible impact on the results. More worrisome is the restriction to ID or 
3D LTE line formation . The atmospheric co nvective motions typically result in C-shaped 



line asy mmetries (e.g . Asplund et al. 2000h . which can thus mimic the presence of 6 Li 



Indeed, Cavrel et al) ( 20071 ) have argued based o n 3D non-LTE calc ulations that the 



Asplund et al.l (|2006( ) results must be reevaluated. Steffen et al.l (|2010f ) have considered 
this in more detail and concluded that the derived 6 Li/ 7 Li ratios in a ID analysis are 
overestimated by typically 0.015, which if true could e x plain the tendency for most stars 
with non-significant detectio ns in the Asplund et al.l (|2006f ) sample to cluster around 

6T ; /7T i~ n m Vnr^nr^nrn iQ + offon of oil Jomnh ovo- 



Li/'Li« 0.01. Furthermore, ISteffen et al.l (|2010l ) argue that by taking this convective 
effect into account, o nly 4 out of 10 stars would remain as > 2cr detections. 

We argue that the ISteffen et al. I (|2010l) results shou ld not be over-inter preted, since 
there are subtle but crucial differences in the analyses. ISteffen et al. I (|2010h rely only on 
the Li line itself to determine all parameters: 7 Li and 6 Li abundances, wavelength shift 
and i ntrinsic line br o adeni ng (rotation in 3D but in ID also macro-/microturbelence) 
while Asplund et al. ( 20061 ) rely on a range of Fe, Ca etc lines to first determine the 
line broadening which is then fixed for the Li profile fitting. Since the other lines are 
also asymmetric, this is accounted for in the profile fitting through a slightly larger 
required broadening, which then compensates for the use of symmetric Li lines in ID. 



Asplund et al.l (|2006f ) also did a 3D LTE analysis for all lines and found very similar 



results as in ID but with much im proved profile fits. W e have also performed 3D non-LTE 
line calculations for Li similar to ISteffen et al.l ( 2010l ) and can in fact exactly reproduce 



their results and their claimed 0.015 effect in 6 Li/ 7 Li when only using the Li line and 
not any other lines. We argue though that this does not make full use of the information 
encoded in the spectra due to the degeneracy between the four fitting parameters. Indeed, 
their procedure is akin to determining the D abundance in QSO absorption system only 
from the Lya line without first resolving the intrinsic velocity structure of the H I clouds 
from other metallic lines. Similarly problematic is to rely on 3D non-LTE for Li but 3D 
LTE for the other lines, in view of the likely pronounced non-LTE effects. The only way 
forward to convincingly demonstrate whether or not 6 Li is present in detectable amounts 
in the atmospheres of halo turn-off stars is to perform full 3D non-LTE calculations for all 
lines considered. We are currently working on this and will report the results elsewhere. 

In summary, it is not yet possible to say that 6 Li has definitely been detected but 
it is definitely too early to say that 6 Li has not been detected. We also note that even 
ISteffen et al. 1 (|2010l) agree that 6 Li is present in HD 84937 and a few other stars. Since 6 Li 
is always destroyed and more so than 7 Li, one would still end up with a very significant 
cosmological 6 Li problem when invoking stellar depletion to solve the cosmological 7 Li 



cosmological Li problem wnen invoking stellar depletion to solve the 
problem (e.g. lAsplund et al.lliool iKorn et ai]l2006t iLind et al.ll2009bT ). 



2.5. Beryllium 

In late-type stars the Be abundance can in practice only be derived from the resonance 
doublet of Ben at 313 nm. In contrast to Li and B, Be is normally present in significant 
amounts in both neutral and singly ionized form due to its high ionization potential, which 
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makes Be less prone to over-ionization. The Be n line formation proceeds under non-LTE 
conditions but because the two dominant non-LTE effects - over-ionization and over- 
excitation - largely c ompensate each other the resulting non-LTE abundance corrections 



are typically minor ([Garcia Lopez et al.1 11995), Both the lower and upper level of the 



doublet are over-populated relative to the LTE prediction. Bel tends to be somewhat 
over-ionized due to a UV radiation excess (J w /B v > 1), which sim ultaneously produce 



an ove r-excitation of the upper level of the Ben 313 nm transitions. iGarcia Perez et al 



(2010) lave investigated the non-LTE line formation of the Ben 313 nm doublet across 



a wide range of stellar parameters. They found that at solar metallicity the two effects 
almost perfectly balance each other with the resulting non-LTE abundance corrections 
amounting to |Alog£B C | ^ 0.05dcx. At low metallicity the non-LTE corrections are 
positive (AlogeBc ~ O.ldex) for T e g ^ 6000K and negative (Alogeec ~ — O.ldex) for 
lower T c g. While the non-LTE effects will not qualitatively change the conclusions for 
the metallicity evolution of Be (Primas and Boesgaard, these proceedings), accounting 
for the new non-LTE calculations will tend to make the slope somewhat shallower. 

The 3D line forma tion of the Be n lines have so far only been investigated for the Sun 
( Asplund et"ai1l2009h . The Be n lines are not very sensitive to 3D effects at least within 



the LTE approximation. The 3D non-LTE case has not yet been investigated but given 
the relatively modest non-LTE effects in ID one would not expect a great difference in 
3D, although this prediction should be verified with detailed calculations. 

2.6. Boron 

As for Be, both over-ionization and over-excitation are at play in the formation of the 
Bi resonance lines at 249.7 nm and 209.0 nm but because Bl is the minority ionization 
stage for ^ 6000 K the e ffects on the resulting line profiles are much more dramatic 



( Kiselman fc Carlsson 19961 ). The over-ionization is largely driven by photo-ionization 



from the excited B I levels since J u /B ll > 1 at the relevant wavelengths. It is also assisted 
by pumping in the B I resonance lines, which produces a sufficient over-population of 
these upper levels. The same pumping results in a substantial increase in the line source 
function S v /B u m J v /B v > 1 for the 249.7 and 209.0 nm lines. Thus, in contrast to the 
case of Be, for the B I resonance lines the two non-LTE effects work in tandem to decrease 
t he non-LTE line strengths su bstantially compared with the LTE prediction. 

Kiselman fc Carlsson (1996) computed non-LTE abundance corrections for the B I lines 



for a grid of MARCS ID model atmospheres. The non-LTE effects grow in size towards 
higher T c g and lower [Fe/H]. Indeed, they found non-LTE abundance corrections as large 
as > +0.5 for typical turn-off halo stars. As both non-LTE processes feed on the UV radi- 
ation field, it is important to include background line opaci ties both for the calcu l ations 
of the photo-ionization rates and the resonance li nes, which Kiselman fc Carlsson! (|1996l ) 



did in a somewhat approximate manner. Recently. iTan et al 



( 2010l ) confirmed the impor 



tance of over-ionization and pumping but they obtained markedly less severe non-LTE 
effects, which can be traced to a more complete treatment of the line-blocking. Their 
estimated non-LTE abundance corrections amount to w +0.3 dex at th e lowest [Fe/Hj , 
when not including the highly uncertain inelastic H collisions according to lDrawin dl968h : 



the corrections would be ~ +0.1 dex when this formula would be app lied without any 
further scaling factor. As shown in Fig. [3l the new non-LTE results o f ITan et all (120101 ) 



make the resulting evolution of B as a function of metallicity somewhat steeper with a 
slope of « 1.3. The reader is referred to Primas (these proceedings) for a discussion of 
the astrophysical implications of su ch a correlation. 

With the exception of the Sun (|Asplund et al.ll2009l) . the 3D line formation of the 



B I lines has not been investigated, neither in LTE nor in non-LTE. In view of the 
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Figure 3. The derived stellar B abundances in halo stars ijTan et al.l 120101 1 as a function of 
[O/H] for LTE (open circles with dashed line the least square fit) and non-LTE (filled circles 
and solid line); upper limits are marked with arro ws. The O abundances have been corrected for 
n on-LTE effects ac cording to lFabbian et al.l (|2009T ) . Also shown are the LT E-based B abundances 
oflTan et alJ (|2010h but corrected according to the non-LTE calculations of Kisclman fc Carlssonl 
l| 19961 ) (open squares and dotted line). Inel astic H c ollisio ns have not been included for the B 
non-LTE results but for O according to the lDrawi 3 l|l968T l formula. 

prominent non-LTE effects and the general expectation that they would become even 
more pronounced in 3D due to the atmospheric inhomogeneities, performing detailed 3D 
non-LTE calculations for B I should have a high priority. 

3. Concluding remarks 

Important progress has been made over the last decade in terms of improving the 
abundance analyses of the light elements in late-type stars by accounting for non-LTE as 
well as 3D effects. The non-LTE line formation of the resonance lines of Li I, Bell, and 
B I share common properties, but the final non-LTE abundance corrections are always 
element specific. Especially differences in ionization potential and the extent to which 
UV radiative transitions contribute to the excitation and ionization balance govern the 
LTE departures. There are however still several outstanding analysis problems. 

Signs are that the new generation of 3D hydrodynamical model atmospheres are in- 
deed highly realistic and thus trustworthy for abundance purposes as they successfully 
reproduce most, if not all, observational tests they have been exposed to. The key de- 
velopment for the future will be to make such 3D models generally available for a wide 
range of stellar parameters and actually to get stellar abundance practitioners to start 
using them routinely. The main obstacle here will likely be to overcome old habits rather 
than say lack of computational power, since it is now possible to perfo rm a 3D-based 
solar /stellar abundance analysis for thousands of lines, at least in LTE (jAsplund et al. 
2009|). On the non-LTE side, the challenges are two- fold: improving the necessary atomic 
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data, especially for collisions, and to couple the non-LTE line formation with 3D model 
atmospheres, which is still very much unchartered territory. This is now certainly fea- 
sible but will require additional (wo)manpower. Only then can we with some degree of 
confidence argue that our analysis efforts match the investments in obtaining impressive 
observational data. As outlined above, such modelling efforts will no doubt lead to many 
reinterpretations of observations with profound implications for our understanding of 
stellar, galactic and cosmic evolution in general and the light elements in particular. 
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